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REPORT  SUMMARY 


During  the  two  years  of  sponsorship  the  theoretical 
understanding  of  Cascade  Charge  Coupled  Device  (C3D) 
electronic  lenses  have  been  greatly  extended.  Specifically 
the  effects  of  the  intermodulation  products,  sidelobe 
resoluion  limitations,  transfer  inefficiencies  are  now  better 
understood  as  a  result  of  this  program.  Additionally  a 
second  generation  device  that  has  both  beam  steer  and  focus 
sections  has  been  fabricated  and  used  to  test  the  validity  of 
the  performance  models  that  have  been  derived.  A  third 
generation  device  has  been  recently  fabricated  that  includes 
three  beam  steer  section  pairs  and  one  focus  section.  This 
advanced  lens  incorporates  the  cumulative  design  rules 
resulting  from  this  research  program.  It  is  anticipated  that 
this  32  channel,  10,000  charge-coupled  element  lens,  will  be 
used  in  a  4.5  MHz  imaging  system  to  be  developed  under  other 
funding. 
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I.  SUMMARY  OP  RESULTS 

This  report  covers  the  work  performed  during  a  two  year 
period  of  sponsorship.  The  research  performed  involved  the 
modelling  and  fabrication  of  an  novel  integrated  circuit  time 
delay  and  signal  summing  element  called  a  Cascade  Charge 
Coupled  Device  (C3D) .  Although  the  device  had  been  invented 
under  previous  sponsorship  of  another  agency  not  enough 
research  had  been  performed  in  order  to  predict  the 
performance  of  a  given  lens  design.  These  two  years  of 
modelling  and  fabrication  sponsorhip  have  enabled  a  much 
greater  understanding  of  the  capability  of  the  C3D  lens  and 
how  it  is  to  be  optimized  for  a  specific  application. 

Although  the  use  of  multi-channel  delay-add  is  generic  to  a 
wide  class  of  electronic  signal  processing  applications,  the 
specific  application  of  interest  to  the  authors  is  the  use  of 
the  device  for  beamsteering  and  focussing  of  signals  for 
acoustic  transducer  arrays. 

Specifically  during  the  first  year  the  emphasis  was  on 
modelling  and  fabrication  process  development  for  the  third 
generation  lens.  The  second  generation  lens  that  was 
anticipated  to  be  fabricated  in  this  first  year  was  actually 
fabricated  prior  to  the  beginning  of  this  program  due  to 
delays  in  the  start  of  sponsorship.  Extensive  modelling  work 
on  transfer  efficiency  was  performed  and  reported  in  the 
previous  semi  annual  reports  of  the  first  year.  In  summary 
it  was  found:  (1)  the  attenuation  effects  of  transfer 
efficiency  are  deterministic  and  insignificant,  (2)  the 
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,  lateral  resolution  and  sidelobe  response  are  to  first  order 

^  not  modified,  and  (3)  for  moderate  values  of  charge-transfer 

inefficiency  the  spatial  shifting  of  the  beamsteering  angle 
is  minor.  Thus  the  C3D  is  largely  insensitive  to  the  major 

*  limitation  commonly  encounterd  (transfer  inefficiency)  in 
charge  coupled  device  applications  today.  A  comprehensive 
coverage  of  this  analysis  is  presented  in  the  Stanford 

*  University  doctoral  thesis  of  John  Shott.  A  complete  reprint 
of  this  dissertation  is  publicly  available  through  the 
University  of  Michigan  Archives  for  a  nominal  reproduction 

*  charge.  Appendix  A  of  this  report  contains  the  analysis  of 
the  fundamentals  of  C3D  performance  found  in  Chapter  5  of 
John  Shott' s  dissertation. 

I  The  fabrication  process  development  of  the  first  year 

concentrated  on  polysilicon  electrode  and  insulating  silicon 
oxide  formation.  Test  fabrication  runs  were  made  to  optimize 
I  further  the  design  parameters  of  channel  length  and 

polysilicon  electrode  conductance.  The  channnel  length 
studies  concerned  experiments  in  electrode  etching 
»  technolgies.  The  previously  used  etchants  were  originally 

developed  for  optical  sensors  with  no  significant  constraints 
on  etching  small  geometries.  The  difficulties  experienced 
I  with  small  geometry  electrode  (5  micron  lines)  etching  were 

associated  with  bubbles  that  formed  on  the  surface  being 
etched.  These  bubbles  radically  enhanced  the  etch  rate  in 
,  their  vicinity  by  over  an  order  of  magnitude.  The  etching 
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process  developed  by  Ernie  Wood  circumvented  this  problem  by 
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scrubbing  the  surface  of  the  wafers  during  etching  with  a 
high  flow  rate  of  nitrogen  gas  bubbles.  The  action  of  the 
bubbles  on  the  surface  increased  the  uniformity  of  the  etch 
by  preventing  bubbles  formed  during  the  etching  process  from 
modifying  one  region  etch  rate  while  regions  without  bubbles 
etched  at  a  different  rate.  The  etching  process  provided 
reproducible  five  micron  line  and  space  etching.  The 
parameters  of  the  etching  process  were  empirically  derived 
using  a  qualitative  success  criteria.  This  fabrication 
research  was  not  performed  under  ARPA  sponsorship. 

Alternative  sponsorhip  was  necessary  as  a  result  of  the 
greatly  reduced  budget  for  the  second  year  of  work  from  that 
which  had  been  originally  planned.  Fabrication  research  was 
supported  by  other  sources  to  make  funds  available  to 
purchase  the  spectrum  analyser  required  for  the  second  years 
work.  The  details  of  these  etching  technologies  will  be  made 
available  in  the  open  literature  in  the  dissertation  of  Mr. 
Wood. 

During  the  second  year  of  funding  additional  modelling 
of  the  inter-modulation  products  was  performed,  the  third 
generation  C3D  lens  design  was  completed  and  fabrication  runs 
were  performed  on  a  device  test  chip.  A  spectrum  analyser 
was  purchased  to  perform  critical  measurements  in  the  the 
verification  of  the  intermodulation  product  model.  This 
instrument  has  proven  invaluable  in  making  spectral 
measurements. 

The  intermodulation  product  modelling  was  concerned  with 
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the  effects  on  performance  resulting  from  the  direct 
connection  of  two  charge  coupled  device  section  operating  at 
different  frequencies  without  the  conventional  lowpass 
antialiasing  filters  commonly  used  in  sampled  data  systems. 
Michael  Eaton  developed  a  simple  analysis  technique  using 
orthogonal  deconvolution  of  the  problem  into  two  parts.  This 
novel  analysis  method  was  described  in  the  third  semiannual 
report  and  in  much  more  detail  in  his  Stanford  University 
doctoral  dissertation.  This  analysis  from  Chapter  IV  of  Dr. 
Eaton's  dissertation  is  repeated  in  in  Appendix  B  of  this 
report.  The  results  of  this  analysis  substantiate  that  the 
abberations  resulting  from  the  intermodulation  process  will 
in  most  practical  cases  not  limit  system  performance. 

The  fabrication  work,  though  extensive,  is  not  as  easily 
summarized  by  few  words.  Perhaps  the  most  impressive  summary 
that  made  is  obtained  by  studying  the  third  generation  C3D 
design  layout  shown  in  Fig.  1  and  the  photomicrograph  of  the 
10,000  element,  32  channel  device  shown  in  Fig.  2.  This  C3D 
represents  the  most  complex  integrated  circuit  device  ever 
fabricated  in  the  Stanford  University  Integrated  Circuits 
Laboratory,  a  laboratory  well  known  throughout  the  world  as 
the  recipient  of  more  outstanding  paper  awards  at  the 
International  Solid  State  Circuits  Conference  than  other 
university  laboratory.  This  complex  C3D  is  more  than 

a  test  structure  to  verify  intermodulation  research  models 
developed  under  the  sponsorship  of  this  grant.  It  is  a 
complete  third  generation  lens  capable  of  imaging  using  4.5 


MHz  ultrasound.  While  this  grant  concerned  the  modeling  of 
intermodulation  distortion  and  other  performance  limits  of 
C3D  lenses,  the  fabrication  techniques,  the  C3D  design 
techniques,  and  the  application  tradeoff  analysis  were 
developed  under  parallel  funding.  The  results  of  this 
research  will  be  made  available  to  other  researchers  in  the 
forthcoming  doctoral  dissertation  of  Ernie  Wood. 

While  funding  of  this  work  has  been  carried  out  under 
the  sponsorship  of  the  Advanced  Projects  Research  Agency,  the 
majority  of  sponsorship  of  the  doctoral  programs  of  Shott, 
Eaton,  and  Wood  were  derived  from  other  sponsor  interested  in 
the  application  of  the  C3D  lens  to  a  specific  application. 
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II.  List  of  Relevant  Dissertations 

A.  John  Shott,  "Charge-Coupled  Devices  for  Use  in 
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IV.  COUPLING  ACTIVITIES  WITH  OTHER  RELATED  RESEARCH 

The  most  significant  coupling  with  other  ARPA  research 
in  related  activies  concerns  our  extensive  interrelationship 
with  the  research  activies  of  Professor  Gordon  Kino  of  the 
Stanford  University  Ginston  Laboratory.  Professor  Kino  has 
long-term  interest  in  the  technologies  used  in  the 
nondestructive  inspection  of  materials  (NDI).  This  research 
effort  and  those  of  Professor  Kino  have  benefitted  through 
both  formal  research  meetings  and  informal  discussions 
between  students  and  researchers. 
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Fig.  2  Photomicrograph  of  32  channel,  10,000  element  C3D  lens 
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APPENDIX  A:  Chapter  5  of  John  Shott's  Ph.D.  Dissertation 
FUNDAMENTAL  LIMITS  OF  C3D  PERFORMANCE 

The  time-delay  accuracy  of  the  cascade  charge -coupled  device  will 
largely  determine  its  ultimate  utility  as  an  electronic  lens  for  an  ul¬ 
trasonic  imaging  system.  ^  '  Time-delay  and 

amplitude  errors  can  reduce  overall  image  quality  by  creating  uncertainty 
as  to  the  position  of  the  focal  point  of  the  lens  or  by  degrading  the 
sldelobe  level  of  the  image.  As  a  result,  this  chapter  examines  the 
aspects  of  the  C3D  electronic  lens  that  will  inherently  limit  the  time- 
delay  and  amplitude  accuracy  generated  by  this  approach. 

The  effects  of  incomplete  charge  transfer  on  the  imaging  system  are 
considered  in  Section  A.  The  remainder  of  this  chapter  discusses  the 
limitations  imposed  by  including  a  number  of  independently  clocked  sec¬ 
tions  on  a  single  device;  also'  described  are  the  effects  of  spurious 
difference-frequency  components  in  the  C3D  output  spectrum  and  the  asyn¬ 
chronous  transfer  of  signal  charge  across  a  frequency  boundary. 

A.  Effects  of  Incomplete  Charge  Transfer 

To  assess  the  effects  of  incomplete  charge  transfer  on  the  perfor¬ 
mance  of  the  20-input  four-section  C3D  electronic  lens  presented  in 
Chapter  IV,  this  device  is  modeled  as  a  4  X  20  array  of  individual  CCD 
delay  lines,  each  characterized  by  the  number  of  bits  in  the  delay-line 
section,  the  clock  frequency  applied  to  the  section,  and  knowledge  of 
Incomplete  charge  transfer  per  transfer  as  a  function  of  clock  frequency. 
Because  the  cascade  charge-coupled  device  (C3D)  lens  is  fabricated  on  a 
single  integrated -circuit  die,  good  matching  of  Incomplete  charge  trans¬ 
fer  vs  clock  frequency  across  the  entire  array  is  expected.  Consequently, 
it  is  assumed  that  a  single  transfer  inefficiency  vs  clock-frequency  re¬ 
lationship  applies  to  all  of  the  time-delay  sections  in  the  C3D  electronic 
lens. 

The  position  of  each  delay  line  within  the  C3D  lens  is  identified 
by  the  subscripts  i,J,  where  i  =■  1,2, 3, 4  refers  to  the  four  clocked 
sections  and  J  ■  -9,-8, . . .  ,-0,-tO,+l, . . .  ,+8,+9  refers  to  the  transducer 


channel  position  numbered  from  left  to  right  across  the  lens  Based  on 

th 

this  notation,  the  time  delay  applied  to  the  j  transducer  channel  (as¬ 
suming  ideal  delay  lines)  is 


AT 


U 


(5.1) 


where 


ducer 


th 

F  is  the  clock  frequency  applied  to  the  i  delay  section  and 

1  th  th 

is  the  number  of  bits  in  the  i  clock  section  of  the  j  trans- 


Recalling  the  dispersion  relations  of  Joyce  and  Bertram  [5.1]  as 
expressed  in  Eqs.  (4.4)  and  (4.5),  attenuation  and  time  delay  AT^ 

for  each  transducer  channel  can  be  defined  as 


1  til 

plied  to  the  i  section  of  delay.  Although  these  expressions  for  at¬ 
tenuation  and  time  delay  are  significant,  the  variation  of  these  quanti¬ 
ties  from  transducer  to  transducer  rather  than  their  absolute  magnitude 


will  determine  overall  image  quality.  As  a  result,  the  following  subsec¬ 
tions  will  examine  the  effects  of  attenuation  and  time-delay  variations 
across  the  array  caused  by  incomplete  charge  transfer.  Because  beam 


*The  subscripts  J  =*  -0  and  j  =  0  refer  to  the  two  center  transducers 
of  an  array  containing  an  even  number  of  transducers.  In  this  context, 
the  "transducer  at  j  +l"  is  located  next  to  the  transducer  at  position 
J  (if  j  »  -0,  then  j  +1  »  0). 


steering  and  focusing  are  generally  independent  of  one  mother,  the 
impact  of  such  a  transfer  on  each  of  these  functions  is  considered  sepa¬ 
rately. 


1.  Beam  Steering 

For  the  beam-steering  section,  Eqs .  (5.2)  and  (5.3)  can  be 
rewritten  as 


f  -  %) 3,3 1  •  \) 4,1  co>  PH] 


(5.4) 


and 


(5.5) 


Maintaining  the  beam-steering  convention  (established  in  Chapter  IV)  that 
the  number  of  bits  in  each  of  the  beam-steering  sections  increases  or  de¬ 
creases  by  an  increment  of  ng  bits  per  transducer  channel,  the  recursive 
relations  that  express  the  attenuation  and  time  delay  of  transducer  chan¬ 
nel  j+1  in  terms  of  channel  j  are 
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and 
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Equation  (5.6)  indicates  that  the  attenuation  at  any  two  adja¬ 
cent  transducers  will  differ  by  a  constant  multiplicative  factor  that 
will  be  dependent  on  the  clock  frequencies  and  on  the  transfer  ineffi¬ 
ciencies  at  those  frequencies  but  will  be  independent  of  the  absolute 
position  along  the  transducer  array.  This  type  of  signal  attenuation  is 
analogous  to  the  tissue  attenuation  discussed  in  Chapter  III.B.l.  Eval¬ 
uation  of  Eq.  (3.43)  for  a  20-element  array  of  1.5-mm  1.5-MHz  transducers 
steered  to  an  angle  of  25°  off-axis  resulted  in  an  element-to-element 
attenuation  factor  based  on  a  tissue  absorption  of  0.898.  If  the  C3D 
lens  is  used  to  steer  the  same  array  to  25°  with  F3  =  5  MHz  (eg =0.002) 

and  F.  »  30  MHz  (e_„  =  0 .025) ,  the  element-to-element  attenuation  co- 
*k  oU 

efficient  is  0.995.  This  demonstrates  that,  even  for  broad  differences 
in  charge-transfer  inefficiency  caused  by  the  large  clock-frequency  vari¬ 
ation  required  to  beam  steer  to  the  extremes  of  the  field  of  view,  the 
change  in  attenuation  across  the  array  resulting  from  incomplete  charge 
transfer  is  substantially  less  than  that  resulting  from  acoustic  path- 
length  differences. 

Comparison  of  Eq.  (5.7)  to  the  beam-steeri:.g  angle/clock-fre¬ 
quency  relation  for  ideal  delay  lines  [Eq.  (4.13)]  reveals  that  the  beam- 
steering  angle  in  the  presence  of  incomplete  charge  transfer  9 ^  is 
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(5.8) 
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where  sin  9  is  the  beam-steering  angle  expected  for  ideal  delay  lines 
and  sin  ££  denotes  the  small  error  caused  by  incomplete  transfer.  Con¬ 
tinuing  with  the  same  clock  frequencies  and  transfer  inefficiencies  as 
before  results  in  9  =  24.62°  (sin  9  =  0.417),  sin  =  -0.0022,  and 

9  »  24.48°. 

€ 

Equation  (5.8)  therefore  indicates  that,  although  incomplete 
charge  transfer  will  alter  the  actual  beam-steering  angle  slightly,  the 
lateral  resolution  and  sidelobe  level  of  the  image  point  will  not  be  af¬ 
fected  by  such  time-delay  variations .  Because  transfer  inefficiency  is 
well  matched  at  any  given  clock  frequency  on  the  single-chip  lens,  the 
device  can  be  operated  at  clock  frequencies  much  higher  than  would  be 
normally  expected,  based  on  the  completeness  of  charge  transfer  at  these 
higher  frequencies.  As  indicated  in  the  left  portion  of  Fig.  5.1,  the 
time  delay  error  due  to  incomplete  charge  transfer  is  directly  propor¬ 
tional  to  the  number  of  bits  clocked  at  each  frequency,  so  good  matching 
of  transfer  efficiency  results  in  a  net  time  delay  error  that  varies 
linearly  across  the  array  for  the  beam  steering  section  shown. 

Implementation  of  electronic  beam  steering  with  a  number  of 
individual  CCDs,  as  shown  in  the  right  portion  of  Fig.  5.1,  will  not  pro¬ 
vide  an  assurance  of  well  matched  transfer  efficiency  between  devices. 

As  a  result,  the  time  delay  errors  developed  across  the  array  due  to  in¬ 
complete  charge  transfer  will  be  randomly  distributed  and  will  tend  to 
degrade  the  sidelobe  suppression  afforded  by  the  imaging  system.  In¬ 
creased  angular  beam-steering  accuracy  and  increased  sidelobe  suppression 
calls  for  CCDs  with  low  transfer  inefficiency  at  high  clock  frequencies 
if  electronic  beam  steering  is  realized  using  individual  CCD  delay  lines. 

In  addition  to  substantially  reducing  the  delay  line  hardware 
required  for  electronically  beam  steered  ultrasonic  imaging,  the  C3D 
single-chip  lens  provides  increased  angular  beam-steering  accuracy  because 
of  improved  matching  of  time  delay  between  transducer  channels .  Further¬ 
more,  improved  matching  of  time  delay  should  make  it  possible  to  achieve 
greater  degrees  of  sidelobe  suppression  with  the  C3D  single-chip  lens  as 
compared  to  Individual  CCD  delay  lines. 
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C30  SINGLE-CHIP  LENS  INDIVIDUAL  CCDs 


Fig-  5.1.  C3D  SINGLE-CHIP  LENS  VS  INDIVIDUAL  CCDs.  (Ideal  Time  Delay 
- ;  Including  Transfer  Inefficiency  - ) 


2.  Focusing 

Based  on  the  above  method,  attenuation  and  time  delay  produced 
by  the  focusing  rather  than  the  beam-steering  section  can  be  written  as 
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In  contrast  to  the  beam-steering  section  where  large  differ¬ 
ences  in  frequency  are  required  to  steer  the  array  to  the  extremes  of 
the  field  of  view,  focusing  over  the  entire  range  requires  clock  fre¬ 
quencies  F^  and  Fg  to  differ  only  slightly.  Because  =  €p2  for 

F,  F0,  it  is  assumed  that  A.  will  also  vary  by  only  a  small  amount 
across  the  transducer  array  and  that  will  nearly  equal  those  val¬ 

ues  expected  for  ideal  delay-line  performance.  To  verify  this  assump¬ 
tion,  ATj  was  calculated  for  the  20-element  array  of  1.5-mm  1.5-MHz 
transducers  using  F^  =  8  MHz  (gp^  «*  0.004)  and  F2  =  10  MHz  (gp2  = 
0.006).  Under  these  conditions,  ideal  delay  lines  would  focus  to  a  range 
of  12.0  cm;  including  the  effects  of  incomplete  charge  transfer  results 
in  a  focal  point  located  at  12.05  cm. 

The  effects  of  incomplete  charge  transfer  on  the  operation  of 
the  C3D  demonstrate  that  the  variation  in  attenuation  across  the  array 
is  deterministic  and  insignificant  when  compared  to  variations  in  atten¬ 
uation  caused  by  tissue  absorption  within  the  body.  Differences  in  time 
delay  across  the  array  resulting  from  incomplete  charge  transfer  shift 
the  focal  point  of  the  imaging  system  in  both  range  and  beam-steering 
angle  but  do  not  degrade  the  lateral  resolution  or  sidelobe  response  of 
the  focal  point.  Even  at  moderately  high  values  of  charge-transfer  in¬ 
efficiency,  however,  spatial  shifting  of  the  focal  point  is  slight  when 
compared  to  lateral  resolution  or  depth  of  focus.  This  insensitivity  to 
incomplete  charge  transfer  is  a  direct  result  of  good  transfer-efficiency 
matching  (at  any  clock  frequency)  across  the  array  which,  in  turn,  is  due 
to  the  C3D  single-chip  lens  approach. 


B.  Effects  of  Asynchronous  Charge  Transfer 

In  the  above  analysis,  it  was  assumed  that  the  C3D  lens  could  be 
modeled  as  a  number  of  independent  delay  lines  and,  therefore,  any  pos¬ 
sible  charge  interactions  at  the  boundary  between  two  sections  driven  at 


» 


different  frequencies  were  neglected.  The  remainder  of  this  chaDter  will 
examine  the  transfer  of  charge  across  a  frequency  boundary.  The  genera¬ 
tion  of  frequency  components  outside  the  desired  signal  passband  and  their 
impact  on  imaging-system  performance  will  be  investigated,  and  considera¬ 
tion  of  the  effects  of  asynchronous  charge  transfer  on  the  signal  will 
include  frequency  shifts  and  time-delay  error. 

** . 

1 .  Generation  of  Difference-Frequency  Components 

Chapter  IV  suggested  that  charge  from  a  CCD  section  being 
clocked  at  frequency  F^  transferred  directly  into  an  adjacent  section 
being  clocked  at  F  results  in  the  generation  of  spurious  signal  com- 
ponents  found  at  the  difference  frequency  |  F^  -  F^  |  and  at  a  number  of 
harmonics  of  this  frequency.  To  simplify  the  analysis  of  asynchronous 
charge  transfer,  consider  Fig.  5.2  where  charge  from  an  electrode  pair 
(1/2-bit  CCD  delay  line)  clocked  at  F.  is  transferred  into  an  electrode 
pair  driven  at  F_.  To  further  simplify,  a  uniform  background  charge 

O 

(rather  than  a  signal  frequency)  of  Q  carriers  will  be  transferred 

o 

into  storage  well  1  each  time  F,  undergoes  a  high-low  (falling  edge) 
transition. 
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Fig.  5.2.  ASYNCHRONOUS  TRANSFER  OF  "FAT 
ZERO"  BACKGROUND  CHARGE. 


As  delineated  by  the  accentuated  clocking  waveforms  in  Fig. 

5.2,  charge  can  be  transferred  from  storage  well  1  into  storage  well  2 

only  when  F  is  high  and  F_  is  low.  As  a  result,  charge  transfer  can 
A  o 


be  initiated  either  at  the  rising  edge  of  F  (when  F  is  low)  or  at 

A  D 

the  falling  edge  of  F  (when  F  is  high).  Under  these  conditions, 
the  time  allowed  for  charge  transfer  can  vary  greatly,  depending  on  the 
timing  of  F  relative  to  that  of  F.  ;  however,  even  if  the  time  al- 
lowed  for  charge  transfer  is  very  short,  the  majority  of  charge  will  be 
transferred  successfully  because  of  the  large  self-induced  fields  pres¬ 
ent  during  the  early  stages  of  the  transfer  [5. 2, 5. 3].  It  will  be 
assumed,  therefore,  that  all  of  the  charge  in  storage  well  1  will  be 
transferred  into  well  2  even  if  the  time  "window"  is  short.  Because  an 
asynchronous  transfer  can  occur  only  three  times  in  any  transducer  chan¬ 
nel  of  the  four-section  C3D  lens ,  the  meaning  of  "all  of  the  charge"  in 
the  above  assumption  can  be  made  less  rigid  because  a  slight  amount  of 
incomplete  charge  transfer  will  not  be  multiplied  repeatedly  as  it  is 
in  the  numerous  "normal"  transfers  encountered  from  input  to  output. 

A  critical  feature  in  Fig.  5.2  is  that,  when  FQ  <  F^,  there 

« 

will  be  intervals  when  a  packet  of  charge  in  well  1  will  not  have  been 
transferred  into  well  2  by  the  time  a  second  packet  of  charge  is  trans¬ 
ferred  into  well  1;  therefore,  the  asynchronous  transfer  of  charge  across 
a  frequency  boundary  results  in  the  nonuniform  distribution  of  the  orig¬ 
inally  uniform  background  charge .  The  repetition  period  (the  inverse  of 
the  fundamental  frequency)  of  this  spurious  signal  component  is  the  low¬ 
est  common  multiple  of  the  periods  of  F  and  F  .  As  a  consequence, 

A  a 

if  F.  and  F„  are  both  multiples  of  F.  ,  the  maximum  repetition 
A  b  me 

period  will  be  1/FJ  which,  in  turn,  implies  that  the  spurious  signal 
inc 

components  will  occur  at  F^  and  its  harmonics.  If  F  and  F_  are 

me  x.  ^ 

also  multiples  of  2F.  ,  the  fundamental  frequency  of  the  spurious  sig- 

ICC 

nal  components  will  be  2Finc  rather  than  F^nc* 

When  >  fa»  a  similar  nonuniform  distribution  of  the  uni¬ 
form  bias  charge  Q  occurs,  which  again  results  in  the  generation  of 
o 

spurious  signal  components  whose  fundamental  frequency  will  be  Finc  i* 

F  and  F  are  multiples  of  F.  .  The  unwanted  signal,  however,  would 
A  B  inc 

be  the  result  of  transferring  an  empty  well  out  of  storage  well  2  before 
it  has  been  filled  by  a  transfer  from  well  1. 

Figure  5.3  is  the  resulting  output  spectrum  of  a  uniform  bias 
charge  that  has  been  transferred  across  a  clock-frequency  boundary  when 
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and  Fg  are  multiples 
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Fig .  5.3.  OUTPUT-FREQUENCY 
SPECTRUM  OF  BIAS  CHARGE  AFTER 
ASYNCHRONOUS  TRANSFER. 


larger  than  the  system  bandwidth 
required  to  achieve  a  given  axial 
resolution  but  sufficiently  small 
to  generate  enough  clock  frequen¬ 
cies  (for  a  given  F  )  to  scan 

max 

the  entire  field  of  view.  The  ultrasonic  center  frequencies  that  fall  in 
the  middle  of  the  regions  free  from  spurious  signals,  thereby  allowing 
for  the  maximum  bandwidth  on  either  side  of  the  center  frequency,  are 
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(5.11) 


Because  the  maximum  bandwidth  under  these  condi- 


B  =  F.  ,  Eq.  (5.11)  provides  the  means  to  select  F,  and 
max  inc  inc 


where  m  s  1,2,3, .. . 
tions  is 

m  (the  number  of  spurious  signal  components  that  occur  at  frequencies 
lower  than  the  carrier  frequency).  Alternatively,  the  ratio  of  maximum 
signal  bandwidth  to  ultrasonic  carrier  frequency  is 
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max  _  inc 

v  =  (m  -  1/2)  F 
o  inc 


(m  -  1/2) 


(5.12) 


To  achieve  a  large  fractional  bandwidth,  therefore,  m  should  be  mini¬ 
mized  which,  referring  to  Eq.  (5.11),  implies  that  Finc  should  be  as 
large  as  possible  for  a  given  center  frequency.  As  observed  in  Chapter 
TV,  Finc  must  be  small  enough  to  enable  the  selection  of  an  appropriate 
number  of  scan  angles  to  image  the  entire  field  of  view.  For  vq  =  1.5 

MHz,  F.  =1  MHz  (m  =  2)  is  appropriate. 

’  inc 


Effects  on  Signal-Frequency  Components — Single  Channel 


The  above  section  examined  the  generation  of  spurious  signal 
components  resulting  from  the  asynchronous  transfer  of  charge  across  a 
frequency  boundary.  For  ultrasonic  imaging  applications,  however,  the 
effect  of  asynchronous  charge  transfer  on  echo  returns  should  be  inves¬ 
tigated  ;  amplitude  accuracy,  time-delay  accuracy,  and  freedom  from  fre¬ 
quency  shifts  must  be  characterized  to  determine  the  ability  of  the  C3D 
lens  to  produce  high-quality  ultrasonic  images . 


Figure  5.4  plots  a  sinusoidal  signal  frequency 


sig 


trans¬ 


ferred  from  a  storage  well  clocked  at  FA  into  a  storage  well  clocked 
at  Fg.  As  indicated  by  the  numbered  arrows,  the  transfer  of  signal- 
charge  packets  across  the  frequency  boundary  occurs  in  three  distinct 
stages . 


Stage  1:  at  the  falling  edge  of  F^ ,  a  packet  of  signal 
charge  is  transferred  into  storage  well  1.  If  a  signal 
charge  is  already  present  in  well  1,  the  new  packet  is 
added  to  it.  Because  there  are  no  storage  sites  preceding 
storage  well  1  as  there  are  in  a  real  CCD,  the  size  of  the 
charge  transferred  into  well  1  will  be  the  size  of  the  sig¬ 
nal  waveform  at  the  time  when  F^  undergoes  a  high-low 
transition. 

Stage  2:  the  accumulated  charge  in  storage  well  1  will  be 
transferred  into  storage  well  2  at  the  rising  edge  of  FA 
(if  Fg  is  already  low)  or  at  the  falling  edge  of  Fg  (if 
Fa  is  already  high) .  This  is  the  actual  asynchronous 
transfer  of  charge.  The  time  allowed  for  this  transfer  is 
variable;  however,  because  most  of  the  charge  is  trans¬ 
ferred  rapidly  as  the  result  of  self-induced  fields,  it 
will  be  assumed  that  asynchronous  transfer  occurs  instan¬ 
taneously.  Any  charge  transferred  during  this  stage  will 
be  added  to  the  charge  already  present  in  storage  well  2. 

Stage  3;  the  accumulated  charge  in  storage  well  2  will  be 
transferred  out  at  the  rising  edges  of  fb-  In  a  real  CCD, 
this  charge  packet  will  continue  to  be  transferred  through 
a  number  of  storage  sites;  in  this  analysis,  it  will  be 
assumed  that  the  packets  leaving  storage  well  2  constitute 
the  final  output  signal. 


A  .  . 


In  Fig.  5.4,  the  signal 

packets  being  transferred  out  of 

well  2  will  depend  not  only  on 

v  ,  F. ,  and  F_  but  also  on 
sig  A  B 

the  timing  of  each,  relative  to 

the  other  two.  The  timing  of  each 

waveform  can  be  defined,  relative 

to  a  common  (but  arbitrary)  time 

reference  (t  in  Fig.  5.4). 
ref 

These  three  waveforms  can  be  com¬ 
pletely  specified ,  therefore ,  by 
a  frequency  and  a  time  delay  (ex¬ 
pressed  as  cycles  of  delay  at  the 
given  frequency)  as 
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Fig.  5.4.  TIMING  DIAGRAM  OF  THE 
ASYNCHRONOUS  TRANSFER  OF  SIG¬ 
NAL  WAVEFORMS. 
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signal:  A  =  A 

o 

clock  1:  C,  »  C 

1  o 

clock  2:  C„  =  C 

i  O 

where 


sgn(x) 


if  x  >  0 
if  x  <  0 


(5.13) 


Because  of  the  asynchronous  character  of  the  transfer  of  charge 
from  storage  well  1  to  well  2,  it  is  difficult  to  derive  an  analytical 
expression  for  the  output  pulse  train  emerging  from  well  2.  For  any  set 
of  frequencies  (v  ,  F  ,  and  F  )  and  relative  timing  information 

S  lg  A  o 

(Agig,  Aa,  and  ^) ,  however,  the  sequence  of  output  pulses  expected 
from  storage  well  2  can  be  generated  (based  on  the  rules  and  assumptions 
presented  earlier  in  this  section)  for  the  three-stage  transfer  process. 
Discrete  Fourier  transforms  (DFTs)  may  then  produce  frequency,  amplitude, 


ft 


and  phase  Information  concerning  the  pulse-output  sequence.  Because  the 
DFT  assumes  that  the  sample  record  is  repeated  ad  infinitum,  the  record 
length  should  contain  an  integer  number  of  cycles  of  the  input  signal  and 
the  two  clocks.  For  a  signal  frequency  of  1.5  MHz  (period  =»  0.667  psec) 
and  clock  frequencies  of  integer  multiples  of  1  MHz  (maximum  periods  1.0 
psec),  for  example,  record  lengths  that  are  multiples  of  2  psec  are  ap¬ 
propriate  . 

Figure  5.5  illustrates  the  output  amplitude  and  phase  of  the 
signal  components  of  a  pulse  train  that  has  undergone  asynchronous  trans¬ 
fer  as  calculated  by  a  4096-point  DFT  [5.4]  of  a  16  psec  record;  the  test 

conditions  were  v  .  s  1.5  MHz,  A  .  =0,  F=  8  MHz,  A  =  0 ,  F  = 10 
sig  sig  A  A  B 

MHz,  and  =*  0.  Apparent  in  the  figure  are  the  clock  difference-fre¬ 
quency  components  described  in  Section  B.l  that  appear  at  2,4,...  MHz, 
the  signal -frequency  component  at  1.5  MHz,  and  the  mixing  products  be¬ 
tween  the  signal  and  various  difference  frequencies  that  occur  at  0.5, 

2. 5, 3. 5,...  MHz.  The  calculated  frequency  spectrum  actually  continues 
beyond  the  10  MHz  limit  in  Fig.  5.5.  The  0  dB  level  was  chosen  as  the 

largest  in  the  output  spectrum— the  F  clock-frequency  component  at  10 

B 

MHZ. 


Fig.  5.5.  OUTPUT-FREQUENCY  SPECTRUM  OF  1.5 
MHz  SIGNAL  AFTER  ASYNCHRONOUS  TRANSFER. 


Fig.  5.6.  PHASE  PLANE  REPRE¬ 
SENTATION  OF  TIME  DELAY  AND 
AMPLITUDE  UNCERTAINTY  CAUSED 
BY  ASYNCHRONOUS  TRANSFER. 


Figure  5.6  presents  the  simulated  amplitude  and  phase  of  the 
1.5  MHz  signal  when  the  timing  of  clock  2  was  shifted  relative  to  that 
of  clock  1  by  varying  Aq  and 
leaving  all  other  parameters 
constant.  The  shaded  portion 
indicates  that  the  phase  of  the 
emerging  1.5  MHz  signal  was  in 
the  range  of  0.3375rt  ±  0.0375jt 
rad  which,  for  a  1.5  MHz  signal 
initially  at  zero  phase,  corre¬ 
sponds  to  a  time  delay  through 
asynchronous  transfer  of  0.1125 ± 

0.0125  nsec.  Simulations  over 
a  wide  range  of  ^  revealed 
that  the  probability  distribu¬ 
tion  of  time  delay  within  this 

region  is  uniform  and  that  no  time  delays  occur  on  either  side  of  the 
specified  time-delay  band.  The  output  amplitude  was  constant  for  all 
values  of  . 

In  the  above  example,  the  mean  time  delay  of  0.1125  psec  is 
the  sum  of  one-half  the  clock  period  of  F^  and  one-half  the  clock  pe¬ 
riod  of  F  •  This  is  the  "ideal"  time  delay  that  could  be  expected, 
based  on  the  observation  that  Fig.  5.4  represents  one-half  bit  of  delay 
at  FA  followed  by  one-half  bit  of  delay  at  Fg.  In  addition,  the  min¬ 
imum  time  delay  through  the  structure  is  equal  to  the  period  of  the  higher 
clock  frequency  (F  =  10  MHz),  and  the  maximum  time  delay  is  equal  to 
that  of  the  lower  clock  frequency.  It  can  be  seen  intuitively  by  exam¬ 
ining  the  waveforms  in  Fig.  5.7  that  the  maximum  time-delay  uncertainty 
is  equal  to  the  difference  between 


the  periods  of  F  and  F 


The 


A  B 

heavy  lines  denote  the  times  al¬ 
lowed  for  charge  transfer  from 
storage  well  1  to  well  2,  and  the 
dashed  lines  indicate  that  wave¬ 
form  2  can  be  translated  in  time 
with  respect  to  clock  1  without 


jsmusnR  or 


Fig.  5.7.  TIME-DELAY  UNCER¬ 
TAINTY  PRODUCED  BY  RELATIVE 
TIME  SHIFTS  OF  CLOCKING  * 
WAVEFORMS . 


affecting  the  positions  of  the  charge  transfers.  The  maximum  amount  of 
shift  (and  thus  timing  uncertainty),  therefore,  is  the  difference  between 
the  two  clock  periods . 

In  the  above  example,  clock  frequencies  were  even  multiples  of 
1  MHz  and  resulted  in  difference  frequencies  spaced  at  2  MHz  intervals. 

If  one  clock  frequency  is  an  odd  multiple  of  1  MHz  and  the  other  is  even, 
however,  the  difference  frequencies  will  occur  at  1  MHz  intervals;  this 
is  important  because  the  1.5  MHz  signal  frequency  can  then  mix  with  the 
difference-frequency  component  at  3  MHz .  One  of  the  mixing  products  will 
occur  at  1.5  MHz  where  it  will  be  added  as  a  vector  to  the  desired  1.5 
MHz  signal,  thereby  resulting  in  greater  time-delay  and  amplitude  errors. 
This  mixing  is  illustrated  in  the  DFT  output  spectrum  in  Fig.  5.8  at 

clock  frequencies  of  8  and  9  MHz.  In 
this  simulation,  signal  frequency  has 
been  offset  to  1.56  MHz  to  emphasize 
that  it  is  indeed  a  mixing  operation 
with  the  3  MHz  difference-frequency 
component  that  produces  the  additional 
signal  component  at  1.44  MHz.  Further 
simulation  reveals  that  the  phase  an¬ 
gle  of  the  mixing  product  will  vary 
throughout  2jr  rad  even  though  that 
of  the  actual  signal  component  may 
only  vary  by 


Fig.  5.8.  OUTPUT  SPECTRUM  OF 
1.56  MHz  SIGNAL,  INCLUDING 
MIXING  WITH  A  3  MHz  COMPO¬ 
NENT. 
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The  time  delay  and  amplitude  uncertainty  of  the  resulting  composite  sig¬ 
nal  (Fig.  5.9)  will  be  far  worse  than  the  time-delay  uncertainty  caused 
by  asynchronous  transfer  in  the  absence  of  such  mixing.  The  addition  of 
the  mixing  product  to  the  desired  signal  creates  uncertainty  in  the  am¬ 
plitude  of  the  composite  signal,  whereas  asynchronous  transfer  only  in¬ 
troduces  phase  uncertainty  in  the  desired  signal. 


It  is  evident  from  the  above  */2 

example  that  the  best  time  delay  and 

amplitude  accuracy  can  be  achieved  if 

the  clock-frequency  increment  F. 

"  me 

is  such  that  none  of  the  spurious 

difference-frequency  components  is 

twice  the  desired  signal  frequency 

v  „  •  With  this  restriction ,  F.  = 
sig  '  inc 

1  MHz  for  a  signal  frequency  of  1.5 
MHz  (see  Section  B.l)  is  no  longer 
appropriate;  however,  when  F^c  -  2 
MHz,  the  mixing  of  signal  and  differ- 
S  ence  frequencies  can  be  avoided . 

Table  5.1  demonstrates  the  results  of  introducing  different 
signal  frequencies  and  phase  angles  into  the  asynchronous-transfer  sim¬ 
ulator  based  on  the  assumption  that  F  and  F  are  accurate  and  sta- 

A  B 

ble  at  all  times.  The  striking  feature  of  this  simulation  is  that  the 
time  delay  of  these  signals  is  a  constant,  independent  of  signal  fre¬ 
quency  and  initial  signal  phase  even  though  it  is  not  equal  to  the  ideal 
time  delay  of  0.1125  jisec. 


Fig.  5.9.  TIME  DELAY  AND 
AMPLITUDE  UNCERTAINTY 
WHEN  MIXING  OCCURS. 


3 .  Effects  on  C3P  Time  Delay  and  Imaging  Performance 

In  the  above  section,  the  discrete  Fourier  transform  was  used 
to  determine  the  effect  of  signal-charge  transfer  across  a  frequency 
boundary.  This  section  will  extend  the  results  obtained  to  estimate  the 
overall  time-delay  accuracy  in  a  complete  C3D  lens  containing  several 
frequency  boundaries  in  a  number  of  parallel  channels.  This  estimation 
will  then  be  applied  to  assess  the  maximum  level  of  sidelobe  suppression 
via  the  CCD  approach. 

The  basic  properties  of  the  asynchronous  charge  transfer  for  a 
single  channel  are  summarized  as  follows. 


There  is  no  shift  in  signal  frequency  whe n  crossing  a  clock- 
frequency  boundary. 


,*i a®..-. 
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•  The  output  amplitude  will  be  a  constant,  Independent  of  the 
f  relative  timing  of  the  two  clock  frequencies. 

'  •  The  time  delay  across  the  frequency  boundary  has  a  maximum 

uncertainty  equal  to  the  difference  between  the  periods  of 
the  two  clock  frequencies.  The  probability  density  function 
within  this  time-delay  window  Is  uniform  for  randomly  se- 
f  lected  values  of  ^  and  £q.  If  FA,  Fg,  ^ ,  and  ^ 

are  exact  and  time  independent,  however,  the  time-delay 
error  will  also  be  time  independent. 


Table  5.1 

TIME  DELAY  ACROSS  A  FREQUENCY  BOUNDARY  VS 
INPUT-SIGNAL  FREQUENCY  AND  PHASE 
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Assuming  that  perfectly  stable  clocking  waveforms  can  be  ap¬ 
plied  to  the  various  sections  of  the  C3D  lens ,  the  time-delay  errors  as¬ 
sociated  with  each  transducer  channel  will  be  identical  and,  as  discussed 
in  Chapter  III,  an  error  common  to  each  transducer  will  not  degrade  the 
quality  of  the  image  point.  A  reasonable  estimate  of  the  degree  of  clock 
stability  necessary  to  ensure  constant  time-delay  errors  is  based  on  the 
fact  that  the  relative  timing  between  two  clock  frequencies  (£T  and  £^) 
cannot  change  by  an  appreciable  fraction  of  the  difference  between  the 
periods  of  F  and  F  during  the  interval  required  for  a  signal  to 
pass  through  the  entire  CCD.  The  ratio  of  the  maximum  time-delay  error 
to  duration  a  signal  is  in  the  CCD  produces  a  good  estimate  of  clock  sta¬ 
bility  S.  If  the  maximum  error  is  assumed  to  be  1  percent  of  the  dif¬ 
ference  between  the  clock  periods ,  then  S  is 


S 


0.01|<1/Fa)  -  (1/F  )  | 

A _ D  , 

(N/2)[(1/Fa)  +  U/Fb>] 


0.02 

N 


Ifa-fbI 
<fa  +  V 


(5.14) 


where  N  is  the  total  number  of  delay  bits.  It  if  *  100,  F  a  10  MHz, 

—5  ^ 

and  PB  -  9  MHz,  for  example,  a  stability  of  S  w  10  is  required, 
which  implies  that 


-  T  tt  ±  io“5) 


FA<t)  =»  Fa(0)(1  ±  10"5)  (5.15) 

Fg(t)  =  FB(0)(1  ±  10“5) 


Such  stringent  demands  on  the  accuracy  and  stability  of  FA  and 
Fq  are  not  easily  met,  particularly  in  a  system  requiring  a  number  of 
such  clocks.  In  addition,  the  rising  edges  of  the  clock  waveforms  in  a 
multichannel  system  must  be  applied  to  all  channels  simultaneously  within 
the  time  limit  of  1  percent  of  the  difference  in  clock  periods.  When' 
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P  a  10  MHz  and  F  »9  MHz,  the  maximum  timing  error  is  o^ly  0.1  nsec; 

A  8 

in  a  parallel-channel  CCD,  the  charging  of  the  wide  storage  electrodes 
can  take  much  longer  if  polycrystalline-silicon  electrodes  are  employed 
in  the  CCD  structure.  For  these  reasons,  it  appears  unlikely  that  the 
signal-frequency  time-delay  error  can  be  held  constant  for  all  elements 
in  the  transducer  array. 

When  clocks  fail  to  meet  the  above  stability  requirements ,  it 
may  be  postulated  that  the  time-delay  error  in  each  of  the  channels  will 
be  the  result  of  a  random  selection  from  the  probability  density  func¬ 
tion  that  describes  the  time-delay  behavior  of  that  particular  frequency 
boundary.  Factors  that  may  cause  a  clock  to  be  considered  unstable  are 

•  insufficient  clock-frequency  accuracy 

•  clock-frequency  "jitter"  caused  by  interactions  between  a 
number  of  different  clock  frequencies 

•  rise-time  variations  resulting  from  the  distributed  RC  na¬ 
ture  of  the  signal -channel  storage  electrodes 

Under  these  conditions,  the  time-delay  error  of  any  signal  channel  at  a 
frequency  boundary  can  be  expressed  as  a  zero-mean  uniform  PDF  random 
variable  whose  width,  expressed  as  a  number  of  cycles  of  signal  frequency, 
is 


slg 


The  standard  deviation  of  the  time-delay  error  (again  as  a  number  of  cy¬ 
cles  of  signal  frequency)  associated  with  the  A-B  frequency  boundary, 
therefore ,  is 


AB 


Vsig  1_  J_ 
■JI5  PA  '  rB 


(5.16) 


In  the  C3D  lens  which  has  several  frequency  boundaries,  the 
time-delay  errors,  associated  with  each  boundary  can  be  considered  inde¬ 
pendent  random  variables ,  based  on  the  same  arguments  that  led  to  the 


I 


conclusion  that  the  errors  associated  with  each  signal  channel  are  inde¬ 
pendent  of  one  another.  Because  the  variance  of  the  sum  of  independent 
variables  is  equal  to  the  sum  of  the  variances  [5.5],  the  total  standard 
deviation  of  time  delay  in  the  C3D  lens  is 


(5.17) 


Using  this  equation  to  determine  the  standard  deviation  of 
time-delay  error  for  any  combination  of  clock  frequencies  coupled  with 
the  tolerance  theory  in  Chapter  III,  the  maximum  level  of  sidelobe  sup¬ 
pression  can  be  predicted  for  each  set  of  clock  frequencies  and,  hence, 
for  each  image  point.  Figure  3.9  has  been  reproduced  in  Fig.  5.10,  with 
the  addition  of  the  range  of 


time-delay  errors  predicted  by 
Eq.  (5.17).  Generally,  image 
points  in  the  center  of  the 
field  of  view  (|8|<12°) 
will  require  medium  or  high 
clock  frequencies  without  much 
frequency  variation.  Most  of 
the  image  points  in  this  region 
should  achieve  >40dB  of  side- 
lobe  suppression  if  proper  apod- 
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ization  is  employed .  In  con¬ 
trast  ,  near  the  extremes  of  the 
field  of  view  (12°  <  |@|<25°) 
where  high  and  low  clock  fre¬ 


Fig.  5.10.  LEVEL  OF  SIDELOBE  SUPPRES¬ 
SION  ACHIEVABLE  BECAUSE  OF  ASYNCHRO¬ 
NOUS  CHARGE-TRANSFER  TIME -DELAY 
ERRORS . 


quencies  are  required  in  the 

beam-steering  section,  the  standard  deviation  of  time-delay  error  tends 
to  be  larger  so  that  maximum  sidelobe  suppression  is  30  to  40  dB  down. 


The  maximum  standard  deviation  can  be  estimated  by  assuming  that  F2  =■ 

F  .  and  that  all  other  clocks  are  running  at  high  frequency  Cat  the 
min 

extreme  of  the  field  of  view).  Equation  (5.17)  can  then  be  approximated 
by 


T.max 


slg 


“s/i"  F 


(5.18) 


min 


When  v  .  a  1.5  MHz  and 
slg 


F  .  =6  MHz ,  this  results  in 

min  ' 


cr_  «  0.1 
T.max 


which  is  a  worst-case  maximum  sidelobe-suppression  level  of  better  than 
30  dB  down  for  a  20-element  transducer  array. 


C.  Summary 

The  fundamental  limits  of  the  imaging  performance  of  the  C3D  elec¬ 
tronic  lens  are  set  by  incomplete  charge  transfer  and  the  asynchronous 
transfer  of  charge  across  clock-frequency  boundaries .  Because  good 
matching  of  transfer  efficiency  across  the  monolithic  C3D  lens  is  ex¬ 
pected,  incomplete  charge-transfer  variations  with  clock  frequency  are 
not  expected  to  affect  adversely  the  resolution  or  sidelobe  level  of  the 
resultant  image  points .  Investigation  of  the  asynchronous  transfer  of 
charge  across  a  frequency  boundary  has  demonstrated  that  the  generation 
of  spurious  signal  components  will  occur  at  the  difference  frequencies 
and  at  their  harmonics  of  the  clock  frequencies;  however,  correct  clock- 
frequency  selection  can  ensure  that  these  spurious  signals  will  not 
interfere  with  the  desired  signal  passband.  Further  analysis  revealed 
that  asynchronous  transfer  introduces  slight  time-delay  errors  in  the 
desired  signal  components.  A  statistical  analysis  of  these  errors  cou¬ 
pled  with  the  tolerance  theory  of  random  arrays,  indicate  that  their 
effect  will  be  small  enough  for  the  C3D  lens  to  produce  images  with  high 
degrees  of  sidelobe  suppression. 
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APPENDIX  B:  Chapter  IV  of  Michael  Eaton's  Ph.D.  Dissertation 

MODELING  AND  ANALYSIS  OF  C3D  ASYNCHRONOUS  CHARGE  TRANSFER 

The  C3D  electronic  lens  is  possibly  the  least  understood  o f  the  four 
major  components  which  comprise  the  preprocessor.  Shott  [iV-l]  has  com¬ 
prehensively  examined  many  aspects  pertaining  to  the  design,  construction 
and  analysis  of  the  device.  However,  analysis  of  asynchronous  charge 
transfer  at  the  internal  clock-section  boundaries  has  not  been  treated. 

The  asynchronous  boundary  is  known  [lV-2]  to  introduce  spurious  intermodu¬ 
lation  products  (IPs)  at  the  sum  and  difference  frequencies.  This  chapter 
fully  analyzes  the  signal  distortions  associated  with  asynchronous  charge 
transfer. 

A.  C3P  Structure 

As  diagrammed  in  Figure  IV- 1,  the  C3D  is 

a  multi-channel  charge-coupled  device.  Each  channel  is  a  multiple  clock 


Fig.  IV-l.  Twenty  channel  Cascade  Charge-Coupled  Device  (C3D) . 
Each  channel  has  four  distinct  clock  sections. 
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section  CCD--the  four  portions  of  the  CCD  are  clocked  at  different  rates. 
The  boundaries  (geometry)  of  the  different  clock  sections  are  tailored  to 
perform  the  beam  steering  and  focusing  required  of  an  electronic  lens. 

In  Figure  IV- 1,  sections  3  and  4  produce  a  linear  time  delay  (beam  steer) 
while  sections  1  and  2  produce  a  quadratic  delay  (focus).  Summation  prior 
to  the  output  is  accomplished  by  removing  the  channel-stop  diffusions  in 
section  1  and  using  a  common  collector  diode. 

Electronic  time  delay  is  performed  by  dividing  an  input  signal  into 
charge  packets  which  are  stored  in  potential  wells  under  the  transfer  elec¬ 
trodes  [lV-3].  The  packets  are  transferred  from  one  electrode  pair  to 
another  much  like  a  digital  shift  register.  The  transferrals  occur  in 
synchrony  with  the  clocking  waveform.  To  relate  the  charge  packets  and 
potential  wells  to  Figure  IV-1,  the  four  clock  sections,  (F^,...,F^)  may 
be  thought  of  as  similarly  shaped  sheets  of  corrugated  tin.  The  troughs 
run  width-wise  across  the  C3D  and  channel  stops  run  lengthwise  to  prevent 
lateral  charge  flow.  C3D  clocking  makes  the  corrugations  ripple  toward 
the  output  causing  charge  transport.  In  each  clock  section  the  rippling 
transport  takes  place  at  a  faster  or  slower  rate  (depending  upon  the  clock 
f  requency) . 

Figure  IV-2  shows  the  two-phase  overlapping-gate  CCD  structure  that 
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Fig.  IV-2.  Charge  transfer  in  a  two-phase  overlapping-gate  CCD. 


is  the  basic  building  block  of  the  C3D  electronic  lens.  Surface-potential 
wells  holding  the  signal  charge  are  diagrammed  for  two  points  on  the  clock¬ 
ing  waveforms.  A  cross-sectional  schematic  of  one  channel  of  the  C3D  ap¬ 
pears  in  Figure  IV-3,  which  clearly  delineates  the  internal  clock  boundaries 
where  asynchronous  charge  transfer  takes  place. 


Fig.  IV-3.  Cross-sectional  schematic  of  the  C3D  lens  showing  only 
one  channel.  Charge  is  asynchronously  transferred  at 
three  Internal  clock  boundaries. 

B.  Circuit  Model  of  the  Asynchronous  Interface 

A  mathematical  circuit  model  of  the  asynchronous  interface  is  an  ex¬ 
tremely  convenient  vehicle  for  developing  intuition  as  well  as  providing 
a  foundation  for  analysis.  The  asynchronous  interface  in  Figure  lV-4a  is 
represented  as  two  abutted  CCDs  with  clock  waveforms  c^(t)  and  c2(t)  . 
Charge  is  transferred  synchronously  inside  either  CCD  section  and  may  be 
represented  as  a  sequence  of  charge  samples  with  1/f^  spacing,  f^  being 

the  respective  clock  frequency.  In  this  case  q(k)  is  the  charge  sample 

th  th  ^ 

in  the  m  storage  well  at  the  k  time  period.  The  time  function  of 

charge  passing  through  the  m*1*  well  is  the  infinite  sequence  {q  (0), 

m 

q^l),  ^(2),  ...}.  Since  the  samples  are  uniformly  spaced  in  time, 
z-transform  calculus  can  be  profitably  applied.  However,  at  the  asynchro¬ 
nous  interface  the  sequence  may  not  be  used  because  of  possible  asyn¬ 
chronous  transfer.  Hence,  the  true  time  function  of  charge  passing  across 
the  interface  must  be  utilized.  An  expanded  view  of  this  interface  is 
shown  in  Figure  IV-4b,  where  clock  waveforms  and  signal  wells  are  defined. 
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Fig.  IV-4.  C3D  asynchronous  interface,  (a)  Abuttment  of  two  CCDs 
resulting  in  asynchronous  interface;  (b)  expanded  view 
of  interface  defining  clock  waveforms  and  potential 
wells  for  circuit  model. 


Each  electrode/well  ("half-bit")  unit  In  this  figure  may  be  analyzed 
using  an  equivalent  circuit  model.  Consider  the  half-bit  to  the  left  of 
the  interface  in  Figure  IV-4b.  The  charge  packets  flowing  into  wel^  may 
be  represented  as  an  impulsive  current  generator,  i^(t),  pulsing  at  a  rate 
f^.  Charge  flowing  out  of  wel^  is  represented  by  a  current  i2(t),  while 
the  charge  stored  in  that  well  may  be  thought  of  as  stored  on  a  capacitor 
C,  (the  MOS  capacitor  above  well^) .  With  these  definitions  the  proposed 
"half-bit"  model  is  shown  in  Figure  IV-5.  The  model  consists  of  the  driv¬ 


ing  function,  i^(t) ,  which  is  the 
synchronous  charge  transfer  into 
well^;  a  controlled  output  cur¬ 
rent  source,  i  (t) ,  and  an  asyn- 
chronous  switch,  S,  which  closes 
whenever  c2(t)  ~  ®^(t)  surpasses 
some  threshold  T. 

A  complete  model  (Figure 
IV-6)  of  the  interface  is  com¬ 
prised  of  two  cascaded  half-bit 
lumps  representing  both  well^  and 


C2(0-C1(t)  >T 
ELSEWHERE 

Fig.  IV-5.  Circuit  model  of  "half¬ 
bit"  clocked  by  c^(t)  in  Fig. 
IV-4b.  Charge  on  capacitor  C 
is  same  as  charge  stored  in 
well  ;  i  (t)  represents  charge 
flowing  Sut  of  we*I^  quanti¬ 
ties  are  defined  in  Fig.  IV-4b. 


well,  of  Figure  IV-4b. 


Input 


» 


Fig.  IV-6.  Two- lump  circuit  model  of  C3D  asynchronous  charge 

transfer.  Switch  S  operates  as  In  Figure  IV-5  while 
second  switch  closes  at  a  rate  f^.  x(t)  is  the  Input 
current;  y(t)  the  output  current. 


x(t)  supplies  samples  at  a  constant  rate  f^.  Switch  S  operates  asyn¬ 
chronously  as  In  Figure  IV-5  with  an  average  rate  of  f  =  min(f^,f2>. 
Current  ix(t)  flowing  through  switch  S  excites  the  second  lump.  The 
switch  in  this  lump  operates  synchronously  at  a  rate  f  resulting  in  a 

Ct 

synchronous  output  y(t) . 

• 

Verification  of  the  accuracy  of  the  model  is  possible  with  the  aid 
of  Figure  IV-7,  which  shows  two  complete  bits  on  each  side  of  the  inter¬ 
face.  This  interface  system  has  only  four  stable  surf ace- potential  con¬ 
figurations,  each  a  combination  of  two  states:  state  0:  c  <f>  >c  <|>0;  and 
state  1:  c  4>„>c  <j>  .  c  $  are  the  four  clock  voltages.  States  11  and 
10  correspond  to  current  source  x(t)  supplying  an  impulse  of  current  to 
capacitor  C  while  switch  S  is  open.  State  01  corresponds  to  switch 
S  being  open  and  x  =  0  with  some  charge  on  capacitor  C.  State  00  is 
the  only  state  possible  where  S  is  closed  and  charge  is  transferred  out 
of  wellj  .  in  this  idealized  model  it  is  important  to  note  that  switch  S 
can  never  be  closed  while  x  is  supplying  current.  The  reason  is  that 
x  ^  0  implies  C^2>C1^1  yi*ure  *v”7)  which  in  turn  implies  c^Cg 
(in  Figure  IV-5).  Therefore  c  -c  s£  0  and  switch  S  is  open  by  virtue 
of  T  being  a  positive  threshold.  In  summary,  the  circuit  model  of  Fig¬ 
ure  IV-6  accurately  represents  all  possible  charge  configurations  at  the 
asynchronous  interface. 


C.  Charge  Transfer  Analysis: 

Orthogonal  Decomposition 

The  system/circuit  designer  In¬ 
terested  in  applying  the  C3D  is 
basically  concerned  with  three  as¬ 
pects:  (1)  frequencies  of  the  spur¬ 
ious  components;  (2)  time  delay  of 
the  signal;  and  (3)  output  amplitude 
of  the  signal.  Thus,  the  difficult 
problem  of  finding  the  general  time- 
dependent  Impulse  response  of  the 
C3D  is  circumvented.  Since  the  de¬ 
sired  signal  component  will  be  fil¬ 
tered  at  the  output,  it  is  suffi¬ 
cient  to  only  determine  the  frequen- 

« 

cles  of  the  interfering  components, 
their  amplitudes  being  unimportant. 
These  interference  frequencies  can 
then  be  avoided  and  the  time  delay 
and  amplitude  of  the  output  signal 
component  are  the  only  remaining 
quantities  of  interest. 

The  prime  stumbling  block  in 
analysis  of  the  model  of  Figure  IV- 6 
is  the  asynchronous  operation  of  switch 
is  defined  as 


C1  c1$2|  c2^1  c2^2 


I 


Fig.  IV-7.  C3D  asynchronous  inter¬ 
face  and  the  four  stable 
charge  configurations.  Con¬ 
figurations  verify  circuit 
model  of  Fig.  IV-5. 


S.  If  a  switching  function,  s(t), 


1,  S  closed 
0,  S  open 

then  the  switch  may  be  considered  as  a  gate  and  s(t)  =  l/2  { 1+ 

3gn[c0( t)-c < t: )  —t] 3 .  The  easiest  way  to  gain  insight  into  the  behavior 
of  s(t)  is  to  consider  a  simple  example  such  as  that  in  Figure  IV-8. 
c.j (t)  is  a  5  Hz.  square  wave,  c^(t)  is  a  4  Hz.  square  wave,  and  T  is 

one  half  of  the  square  wave  amplitude.  In  this  case  the  gating  function 

16 


c2<‘>:  j~Lrijnj-LJij'uxrm_ri^ 


ORTHOGONAL 
DECOMPOSITION 


INPUT:  < 

TRANSFER: 

OUTPUT: 


a2  0  *3  *4  *5  *6  0 


Fig.  IV-8.  Decomposition  of  gating  function  s(t)  into  orthogonal 
basis  vectors.  c2(t)  is  a  5  Hz.  square  wave,  c^(t)  a 
4  Hz.  square  wave?  s(t)  has  a  1  Hz.  period  which  is 
the  greatest  common  denominator  (GCD)  of  4  and  5. 


s(t)  is  given  by  c2(t).c^(t),  the  dot  being  Interpreted  either  as  multi¬ 
plication  or  logical  ANDing.  The  overbar  means  the  logical  complement. 

For  ease  of  graphing,  c  and  cT  are  plotted,  and  s  =  1  iff  both  c  =  1 

2  i  t  2 

and  c^  s  1.  Assuming  instantaneous  charge  transfer  [lV-4],  the  leading 
edges  of  s(t)  mark  the  points  when  the  asynchronous  transfer  occurs. 

Input  pulses  occur  on  the  falling  edges  of  c^(t)  while  output  pulses  are 
on  falling  edges  of  c2(t).  The  figure  also  shows  a  sequence  of  periodic 
inputs,  the  transfer  sequence  across  the  Interface,  and  the  periodic  out¬ 
put  pulses  resulting  from  the  input  sequence. 

Analysis  is  facilitated  by  noting  that  s(t)  may  be  decomposed  into 
four  orthogonal  components,  each  periodic  with  a  common  frequency.  This 
is  given  by  the  greatest  common  denominator  (GCD)  of  4  and  5,  which  is  1  Hz 
Since  all  components  have  the  same  period,  the  sum  is  periodic  with  rate 
GCD  also.  Thus,  the  asynchronous  interface  is  a  time-varying  system  with 
periodic  parameters. 
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4 


4 


3 


4 


I 


* 

i 
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The  general  analysis  of  time-varying  systems  is  unwieldy  but  since 
the  asynchronous  interface  is  a  linear  system  with  periodic  parameters, 
we  may  employ  the  simplifying  theory  of  Floquet  [lV-5].  Briefly,  Floquet 
theory  states  that  if  an  input  of  the  form  exp(j2jtfQt)  is  applied  to  a 
linear,  time-varying  system  with  periodic  parameters  of  frequency  GCD,  the 
output  has  the  form 

00 

2  Bk  exp[  J2  jt(  f  Q+kGCD)  t  ] 

This  situation  is  shown  in  Figure  IV-9a.  In  the  case  of  the  C3D  asynchro¬ 
nous  interface,  a  sinusoidal  input  is  sampled  at  frequency  f ^  (the  fre¬ 
quency  of  clock  c^(t)).  The  corresponding  input  and  output  spectra  are 
schematically  diagrammed  in  Figure  IV-9b.  Assuming  both  a  real  input  sig¬ 
nal,  cos(2)tf  t),  and  dc  bias  ("fat  zero")  charge,  the  output  spectrum 
o 

contains  components  at  frequencies  of  the  form 


±  (Zf  +mf  +nGCD) 
o  l 


where  m  and  n  are 

integers  and 

Z  =  one  or  : 

zero. 

...  t 

t 

t 

t  V  ..  , 

f0 

freq. 

fo'f» 

*0 

f0+'1 

f0-2f,  freq. 

...  f 

t 

t  t  —  . 

t  t 

t  — 

VGCD 

f0 

f0*GCDf0*2GCD 

f0  VGCD 

v*« 

freq. 


(a) 


(b) 


Continuous  Input  Sampled  Input 


Fig.  IV-9.  Input  and  output  of  linear,  time-varying  system  with 
periodic  parameters  for  (a)  continuous  input  and  (b) 
sampled  input.  GCD  a  parameter  rate;  fQ  =  signal 
frequency;  f  =  sampling  frequency. 
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If  f ^  s  pGCD  where  p  is  an  integer,  the  interference  frequency  is 

±[Zf  +(mp+n)GCD].  Interference  with  the  desired  signal  component  occurs 
o 

when 

f  =  -  C  f  + ( mp+n) GCD] 
o  o 

OP 


f  =  *  (mp+n) GCD 
o 


which  has  solutions  kGCD/2.  Therefore  the  objective  of  no  interference 
translates  into  the  requirement 


f 

o 


k  any  integer 


For  example,  if  'GCD  =  1MHz  then  f  A  0.5,  1.0,  1.5,  2.0,  2.5,  ...  MHz. 

o 

Thus .  GCD  =  1MHz  Is  a  reasonable  choice  for  f  =  2.25  MHz ,  a  standard 

.  o 

ultrasound  frequency.  It  is  clear  from  the  above  that  the  maximum  band¬ 
width,  BW,  of  the  ultrasound  burst  is  restricted  by  BW  <  GCD/2  from  which 
the  following  design  inequalities  may  be  obtained: 


GCD  >  2BW 

k  W  i  <fi,  V  -  <f.  -  f  •  fo  +  f  >  k  any  integer 


For  example,  in  modern  pulse-echo  ultrasound 


BW 


f 

o 

2 


then 


r  GCD  >  f 
GCD  >  f 

o 


+ 


BW 

2 


so  the  value  of  GCD  which  simultaneously  satisfies  both  inequalities  is 

GCD  >  4  f  .  To  ensure  some  guard  band,  BW  —  f  would  actually  be  used 
4o  o 


with 


GCD  >  |  f 


o 


for  high  resolution  systems. 


Analytic  analysis  of  the  amplitude  and  time  delay  of  the  desired  sig¬ 
nal  component  is  possible  through  use  of  orthogonal  decomposition  of  s(t) 
Since  the  interface  lodel  of  Figure  IV-6  consists  of  two  identical  circuit 
lumps,  each  may  first  be  analyzed  independently  then  later  combined.  The 
analysis  proceeds  by  considering  the  exciting  source  to  be  continuous  in¬ 
stead  of  pulsatile.  By  virtue  of  linearity,  the  sampled  nature  of  the  in¬ 
put  is  constructed  as  a  superposition  of  sinusoidal  Inputs. 


Consider  only  the  first  circuit  lump  of  Figure  IV-6.  Through  appli¬ 
cation  of  the  orthogonal  decomposition  of  s(t),  1^(0  is  a  train  of  im¬ 
pulses,  which,  for  our  purposes  may  be  represented  as 


ix(t) 


L 

2  a  (t) 
1=1  1 


6<t  “  g5d  ~  ti) 


L 

=  2  a  (t)  GCD  comb  [GCD(t-t  )] 

i=l 


where  the  L  basis  vectors  have  leading  edges  at  the  times  .  Suppose 

the  continous  input  x(t)  =  exp  (J2«f  t)  is  applied.  If  (t.  }  is  the 

O  K 

sequence  of  switch  closures  then  a^(tk)  equals  the  strength  of  the  Impulse 

at  t=tj£  which  also  equals  the  Integral  of  input  current  between  tfc  and 

t.  , .  There  are  only  i  different  Intervals  between  switch  closures,  call 
k-1 

these  Aj .  Then 


w 


x(t)dt  = 


exp( j2rtfQt)dt 


*  J2 j  exP<J2«f0tk>  -  exp[j2«fo(tk-Ai)] 


=  At  exp  <-jfffQA1)  sinc(foAt)  expU2>rfotk) 


1  <t)  *  ?  GCDA.  exp(-J«f  A.)ainc(f  A.)  exp(J2?tf  t)  comb  [GCD(t-t  )] 
x  11  OlOl  O  1 


Thus 


t 


and  the  Fourier  transform  of  i  is 

x 

exp(-j2itt  f) 

III,)  .  I00D41«P(-jrfo41).l»=(fol1H(f-fo).  - sar—  ««*  (^) 

(IV-1) 

Equation  IV-1  is  the  general  formula  for  the  spectrum  of  i^tt)  when  the 
input  is  a  continuous  sinusoid.  However,  in  the  real  device  x(t)  is  a 
sinusoid  sampled  at  rate  f^,,  i.e. 


£  exp[ J2 it ( f o+nf  j)  t  ] 


(IV-2) 


Plugging  (IV-2)  directly  into  (IV-1)  yields 


5{ix}  *  2  GCDAt!  2  exp[-j«(fo+nf1)dl]  sinc[(fo+nf 1)Ai]6(f-fQ-nf j) 
1=1  (  n=»  -*» 


exp  (-J2xt  f) 

*  - rpr: -  comb(— — ) 

GCD  GOD 


(IV-3) 


Now  suppose  f^  s  pGCD,  p  an  Integer,  then  the  transform  of  IV-3  evaluated 

at  f  s  f  is 
o 


jtg 


£  GCDAt  £  exp{  jjt[nf  1(2t1-^l)-fo^l]}sincC(fo+nf  (IV-4) 


f=f  1=1  n=-« 

o 


Analysis  of  the  second  circuit  lump  is  exactly  analogous.  The  output 
signal  has  the  form 


y(t)  =  £  b(t)5(t-  —)«  b(t)  f2  comb(f  t) 

-«  2 


If  ix(t)  is  the  continuous  sinusoid 


ix(t)  a  exp(J2«fot) 


1 


t 


f 


f 


then 


f 


f{y}  =  exp(-Jj sine  (^B)6(f-fQ)* 

2 


1  u  x 
7-  comb  (— ) 

2  z2 


The  actual  input  to  the  second  circuit  lump  is 

L  » 

i  (t)  =  2  2  exp(-J2fft.nCCD)exp[j2jt(fo+aGCD)t] 

i=l  n=  -<=° 


♦ 


t 


♦ 


and  the  output  spectrum  at  f 


f  is 
o 


L  ® 

2  2 

1=1  n=  -® 


V'S 


exp[  j«(2tinf2-fo/f2]sinc(~) 

2 


(IV-5) 


Therefore  the  frequency  response  at  the  signal  frequency 
product 


H(fo) 


•f(y] 


f 

o 


is  the 


The  conclusion  is  that  the  analytical  result  is  mathematically  tractable 
but  is  too  Involved  for  design  purposes.  A  general  observation  is  that 
the  sine  terms  indicate  that  the  response  falls  off  as  the  signal  frequency 
increases— the  asynchronous  Interface  behaves  as  a  lowpass  filter. 

A  heuristic  approach  will  be  taken  to  obtain  results  which  are  readily 
useful  for  design  purposes.  Time  delay  through  the  asynchronous  interface 
is  investigated  first.  Shott  [lV-6]  has  reported  the  ideal  time  delay 
through  the  interface  to  be  +  y-) •  He  also  found  that  there  is  a 

time  delay  Jitter  which  depends  upon  ?he  relative  phasing  of  the  two  clocks, 
c^t)  and  c2(t).  The  Jitter  can  be  seen  in  the  example  of  Figure  IV- 10. 

In  this  example  f2  =  4 Hz  and  f^  =  1  Hz.  The  ideal  delay  equals  0.625  sec. 
In  Figure  IV- 10a  the  relative  clock  phasing  has  been  adjusted  to  produce 
a  time  delay,  rd  =  0.75  sec.  Two  asynchronous  transfers  per  l/GCD  are  evi¬ 
dent,  however  only  the  first  is  an  active  transfer.  If  it  is  attempted  to 
Increase  the  delay  further  by  slightly  advancing  the  phase  of  c^  as  in  Fig¬ 
ure  IV-lOb,  the  delay  instead  Jumps  to  a  minimum  value  of  =  0.5+g  sec. 
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a0  a1 

♦  ♦ 

♦  * 


♦  * 
♦  ♦  ♦  ♦ 


INPUT: 

TRANSFER: 


**************  OUTPUT: 
0  0  «00  0  0  a,  0  0  0  a20  0  0 


♦  ♦  ♦  ♦  ♦  ♦  *  *  * 


0  aQ  0  0  0  a,  0  0  0  a2  0  0  0  a3 


j  0.75*  |— 


(a) 


— U-  0.5** 
(b) 


Fig.  IV-10.  Relative  phasing  of  clocks  c  and  c alters  the  signal 
delay  through  the  asynchronous  Interface,  (a)  Delay  = 
0.75  sec;  (b)  delay  =  0.5+e  sec. 


This  happens  because  a  third  asynchronous  transfer  appears  at  an  earlier 
time  than  the  original  transfer,  expediting  the  throughput  of  the  charge 
sample.  The  third  transfer  la  a  result  of  the  large  difference  between 
f  and  f_.  Careful  scrutiny  of  the  effects  of  changing  the  relative 

phase  of  c.  and  c  In  Figure  IV" 10  lead  to  the  conclusion  that  the 

J»  z 

peak-to-peak  Jitter  is  equal  to  the  difference  between  the  clock  periods 
reduced  by  a  factor  proportional  to  the  relative  difference  between  f ^ 
and  f  .  Specifically,  If 

fL  g  min^,  f2> 
fH  4  -«(f1,  f2) 

T.(nomlnal)  g  5(7-  +  7-) 
d  2  fx  f2 


then 


1,1  1  .  ±  GCD  1.1  1  . 

Td  *  2^f  +  f  ^  t  -f  ~2^t  "  f 
L  H  H  L  L  H 


-  2(i_  +  L.)  ±  gcd 

“  2'f  t  2ft 
L  H  H  L 


or 


Td  =  Td<nomlnal) 


GCD  - 

w 


(IV-6) 


Thus  the  Jitter  Is  usually  substantially  smaller  than  that  reported  by 
Shott  [lV-6]. 

Since  the  above  formulation  was  developed  from  a  heuristic  standpoint, 
a  computer  simulator  was  written  to  test  for  correctness.  The  asynchro¬ 
nous  transfer  algorithm  Is  from  Shott  [lV-7]  and  is  depicted  in  Pascal¬ 
like  notation  below. 


(*  tfc  are  the  switch  closure  times  *) 

well  =  well  =  0  (*  T  =t  Input  period;  T_  =  output  period  *) 
X  2  X  2 

m  a  n  *  k  a  0 

A  =  min  (TlfT2)/100  (*  a  "small"  initial  phase  *) 

repeat 

time  s  mln(tk,  +  A,  nTg) 

If  time  a  mT^  +  A  then 

wellj^  —  wellj  +■  Input  (m) 
m  —  succ(m) 

If  time  =  tfc  then 

well2  —wellg  +  well^ 
well^  —  0 
k  —  succ(k) 

If  time  =  nTg  then 

output(n)  — wellg 
weU2-0 
n  a  succ(n) 


until  forever 


The  results  of  the  simulator  on  a  wide  variety  of  clock  frequencies, 


relative  phasing  and  signal  frequencies  Indicates  that  equation  IV-6 
accurately  predicts  the  bounds  on  the  peak-to-peak  jitter.  If  the  relative 
clock  phase  has  a  uniform  distribution  then  the  jitter  is  also  uniform  on 

(-  gv-:  ",yj- 1-) .  Random  time  delay  jitter  may  be  eliminated  entirely  by 
phase  locking  the  C3D  clocks. 

The  final  topic  concerning  asynchronous  Interface  is  the  variation 

of  the  output  signal  amplitude  versus  both  clock  and  signal  frequencies. 

The  equation  s  =  e. .c  implies  that  there  is  at  least  one  asynchronous 
1  2 

transfer  every  l/f^  seconds.  Because  of  this,  the  half-bit  (circuit 

lump)  which  clocks  at  f  will  operate  essentially  synchronously,  that 

is,  charge  packets  are  never  added  together  before  being  clocked  out.  The 

signal  therefore  undergoes  pure  time  delay  in  the  half-bit  section  which 

is  clocked  at  f  _  giving  rise  to  two  different  cases : 

L 


ix(»> 


T~1  i 

lii 


Case  (i)  f^  a  fj 

Referring  to  Figure 
IV-11,  the  charge  samples 
ix(t)  are  essentially  de¬ 
layed  versions  of  x(t) , 
so  for  amplitude  analysis 
the  first  lump  may  be 
dropped  and  ix(t)  becomes 
the  input  signal.  The  only 
portion  of  the  Interface 

model  influencing  the  output  signal  amplitude  is  the  second  lump,  shown 
dotted  is  Figure  IV-11.  This  circuit  is  equivalent  to  a  zero-order  hold 
(ZQH)  device  with  a  ’’hold"  time  of  l/f2,  and  frequency  response  [lV-8] 


Fig.  IV-11.  Output  amplitude  analysis. 

fl  *  fL  -  Bln(f i  >f2^  ’  7116  seCond  clr“ 

cult  lump  is  equivalent  to  a  zero- 

order  hold  (ZOH)  circuit. 


|H(fo>|  »  sine  (-2) 
2 


(IV-7) 


The  two  operations  of  sampling  the  model's  capacitor  voltage  and  evalua¬ 
ting  the  sampled  spectrum  at  f  *  f  effectively  "undo"  one  another 
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f 


f 


t 


t 


♦ 
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* 
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with  the  result  that  equation  IV-7  is  the  overall  frequency  response  of 
the  asynchronous  interface  when  f  =  f  . 

Li  X 

Case  (ii)  fL  =  f2 

In  this  case  input  pulses  arriving  at  the  faster  rate  f  (=f  )  are 

H  I. 

averaged  in  both  of  the  capacitors.  However,  those  not  averaged  in  well^ 
are  averaged  in  wellg  and  vice  versa  with  the  net  effect  that  the  input 
samples  are  averaged  over  a  period  of  l/f2  and  tIlen  dumped.  Once  again 
this  constitutes  a  ZOH  device  with  1 /f2  "hold"  time  and  resulting  response 
|  H(f  Q)  |  =  sinc(f  7116  overa11  amplitude  frequency  response  is  the 

same  in  both  cases  and  equation  iv-7  adequately  describes  the  response  of 
the  asynchronous  Interface  for  all  combinations  of  clock  and  signal  fre¬ 
quencies  as  well  as  relative  clock  phase. 

The  computer  simulator  was  again  used  to  check  the  accuracy  of  equa¬ 
tion  IV-7  and  the  results  are  plotted  in  Figure  IV- 12.  Very  good  agree¬ 
ment  was  obtained  in  all  cases;  the  greatest  deviation  from  sinc(f  /f  ) 

O  6 

behavior  was  obtained  using  small  values  of  f^GCD  and  f2/GCD  where 
there  are  fewer  input  samples  averaged  over  the  l/GCD  fundamental  period. 

D.  Summary 

The  important  design  formulas  developed  in  this  chapter  are  presented 
together  here  for  easy  reference. 

Undesirable  signal  frequencies: 

f  si  k  —5-  ,  k  any  integer 

O  J 

Restriction  on  the  greatest  common  denominator  (GCD)  of  clock  frequencies: 


GCD  >  2(BW) 


k 


GCD 

2 


i 


<frV  =  (1o 


BW 

2 


i 


Pig.  IV- 12. 


Amplitude  frequency  response  of  C3D  asynchronous 
interface  for  various  clock  frequency  combinations 
(fj,f2).  Solid  lines  are  sine  (f  Jt 2>  ;  plotted 
points  are  output  from  discrete-transfer  computer 
simulator  of  Interface. 


Nominal  time  delay  of  asynchronous  interface: 


Td(nominal)  =  +  y-) 

X  2 


Time  delay  jitter  using  unlocked  clocks: 


t j  =  Tj (nominal)  [l  *  v00?  ] 


Frequency  response  of  asynchronous  Interface: 


|H<V|=sinc  (f5) 

2 

Several  observations  may  be  drawn  from  these  formulas.  First,  C3D 
clocking  waveforms  should  be  phased- locked  to  prevent  time  delay  errors 
which  were  shown  in  Chapter  III  to  degrade  sldelobe  performance.  If 
phase- locking  is  Impractical,  then  high  clock  frequencies  are  preferred 
which  are  noncomposite  numbers,  i.e.,  their  GCD  is  small.  However,  GCD 
cannot  be  decreased  too  much  or  the  ultrasound  bandwidth  will  be  restricted— 
as  usual,  an  engineering  compromise  must  be  struck.  Finally,  cascaded 
multiple  sections  with  asynchronous  boundaries  must  respect  the  penalty 
paid  in  SNR  as  the  number  of  sections  increases.  In  this  case,  equation 
XV- 7  is  useful  for  C3D  design  purposes  since  each  boundary  behaves  as  an 
independent  transfer  function. 
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